Skyrmions without Sigma Models in Quantum Hall Ferromagnets 



A.H. MacDonald\ H.A. Fertig^, and Luis Brey^ 
^ Department of Physics, Indiana University, Bloomington IN 4^405 
^ Department of Physics and Astronomy, University of Kentucky, Lexington KY 40506-0055 
^ Instituto de Ciencia de Materiales (CSIC), Universidad Autonoma C-12, 28049 Madrid, Spain 

(February 1, 2008) 



On 
OS 



O 

O 



> 
o 

oo 
O 
O 



o 



We report on a microscopic theory of the Skyrmion states 
which occur in the quantum Hall regime. The theory is based 
on the identification of Skyrmion states in this system with 
zero-energy eigenstates of a hard-core model Hamiltonian. 
We find that for A^^ orbital states in a Landau level, a set 
of Skyrmions states with orbital degeneracy A^^ — K and spin 
quantum number S = N/2 — K exists for each nonnegative 
integer K. The energetic ordering of states with different K 
depends on the interaction potential. 
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The ground state of a two-dimensional electron sys- 
tem (2DES) in the strong magnetic field regime of the 
quantum Hall effect is ferromagnetic at certain values of 
the Landau level filling factor v. {v = N /N^ where N 
is the number of particles, = AeB/hc = A/{2tt£'^) 
is the orbital degeneracy of the Landau levels and £ is 
used below as the unit of length.) The simplest example 
of a quantum Hall ferromagnet (QHF) occurs at u = 1 
where the ground state is a strong ferromagnet with total 
spin quantum number S = N/2. Phenomena associated 
with spontaneous magnetic order are open to experimen- 
tal study in QHF's, despite the strong magnetic fields, 
because the Zeeman coupling to the electronic spins is 
smaller than other typical energy scales of the system and 
can even be tuned to zero, for example by application of 
hydrostatic pressure to the host semiconductor. QHF's 
have the unusual property, first identified in finite-size ex- 
act diagonalization studies and dramatically evident 
in recent Knight shift spin-polarization measurements , 
that S can be sharply reduced |^ , in appropriate circum- 
stances even to S" = 0, by the addition or removal of a 
single electron. This behavior may be explained Q using 
a non-linear cr model [NLa) continuum field theory de- 
scription of QHF's. In two dimensions, the NLa model 
supports spin texture excitations, known as Skyrmions, 
that carry a unit quantized topological charge |q] . It has 
been conjectured and under appropriate circum- 

stances it may be explicitly proven Q, that for QHF's 
the electrical charge of a Skyrmion is equal to the prod- 
uct of the ground state filling factor and the topological 
charge. This implies that Skyrmions will be present in 
the ground state for v close but not equal to 1, explaining 
P| the reduction in the spatially averaged moment. 



While the NLa model provides a pleasing qualitative 
explanation of the spin-polarization experiments, it can- 
not be used to address quantitative issues. It is valid 
only for slowly varying spin-textures, while the Zeeman 
coupling at experimentally relevant fields favors small 
Skyrmion states with only a few reversed spins and rel- 
atively rapid variation in the spin-moment orientation. 
Recent microscopic Hartree-Fock estimates of the 
optimal skyrmion size agree well with experiment 
further bolstering evidence for unit charge, large spin, 
quasiparticles of the v = \ ferromagnetic ground state. 
In this Letter we report on an independent, fully micro- 
scopic, picture of QHF Skyrmions. In addition to giving 
an alternate physical picture of these exotic quasiparti- 
cles, our approach has the advantage that we are able to 
precisely determine the quantum numbers and multiplic- 
ities of all Skyrmion states. 

The approach we have taken is in the same spirit as 
the illuminating outlook on the spin-polarized fractional 
quantum Hall effect which arises from appropriate hard- 
core model Hamiltonians ||J^. As discussed for the case 
of interest below, these models have zero energy many- 
particle eigenstates which are often known analytically, 
are separated from other many-particle states by a fi- 
nite gap, and have a degeneracy which increases with 
increasing N . An incompressible state responsible pC| ] 
for a quantum Hall effect transport anomoly in such a 
model is the nondegenerate maximum N zero energy 
eigenstate. The zero energy eigenstates at lower densi- 
ties constitute the portion of the spectrum which involves 
only the degrees of freedom of the, in general fractionally 
charged quasiholes of the incompressible state. It is 
assummed that the difference between the model Hamil- 
tonian and the true Hamiltonian is a sufficiently weak 
perturbation that the quasihole states are still well sep- 
arated from other states in the Hilbert space, although 
accidental degeneracies will be lifted in the spectrum of 
the true Hamiltonian. Here we apply this approach at 
V = 1. Our principle results may be summarized as fol- 
lows. The zero energy N-fermion eigenstates for a sin- 
gle hole in a Landau level may be mapped to a set of 
N-boson states in which the bosons are allowed to oc- 
cupy only four single-particle states. Single-hole states 
exist with total spin number S = N/2 — K for each non- 
negative integer K and in the absence of disorder and 
Zeeman coupling have degeneracy g — gorb dspin where 
gspin = 25* + 1 is the spin-multiplicity, and the orbital 
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degeneracy gorb = N + 1 — K. 

For our analysis we use the symmetric gauge in which 
the single-particle orbitals in the lowest Landau level 
are 

(^;;rTw)V?exp(-|z|V4), (1) 

where |l2| m = 0, 1, • • • , — 1 , z = x + iy, and x and 
y are the Cartesian components of the two-dimensional 
coordinate. We study here the hard-core model for which 
the interaction is ||l3|] : 

V = A7rVo Y^ S^^Hf^-r,) (2) 

i<j 

At strong magnetic fields the low-energy Hamiltonian is 
simply the projection of this interaction onto the the low- 
est Landau level |^ . Many-particle wavefunctions which 
are zero energy eigenstates of this Hamiltonian must van- 
ish when any two-particles are at the same position and 
must therefore have the difference coordinate for each 
pair of particles as a factor: 

nz,x]^[l[{z^-ZJ)]^B[z,x]■ (3) 

i<j 

We note that the each complex coordinate appears to the 
power A^— 1 in the factor in square brackets in Eq. (||) and 
that this factor is completely antisymmetric. It follows 
that ^'b[z] must be a wavefunction for N bosons and 
that these bosons can be in states with angular momenta 
from to N^f, — N. This simple observation leads to the 
conclusions we reach below. 

First we consider the case of a filled Landau level, 
N = N^. In this case all bosons must be in orbitals 
with TO = 0. ^'b[2;,x] must then be proportional to a 
symmetric many-particle spinor and therefore have total 
spin quantum number S = N/2. The orbital part of the 
fermion wavefunction can be recognized as the Slater de- 
terminant with all orbitals from to = to m = iV^ — 1 
occupied. We are able to conclude that the ground state 
is a strong ferromagnet with no orbital degeneracy. The 
ease with which this conclusion can be reached contrasts 
markedly with the case of the Hubbard model where 
enormous effort has yielded relatively few firm results 
1^^ . When Zeeman coupling is included in the Hamil- 
tonian the ground state will be the member of this mul- 
tiplet for which all spins are aligned with the magnetic 
field, i.e. the state with ^ S ^ N/2. 

Our primary interest here is in the elementary charged 
excitations of the ferromagnetic — 1 ground state which 
occur at = A^^ ± 1. For the case of a single hole 
{N = N^ — 1) the lowest energy states are the zero-energy 
eigenstates of the hard-core model. We will see that they 
are the quantized version of the charged Skyrmion states 
in the NLa model description of quantum Hall ferromag- 
nets. From Eq. (^) it follows that these states can be 



mapped to spin-1/2 A' boson states where the bosons 
can have angular momentum equal to or 1. In un- 
derstanding these states it is helpful to start from the 
state with boson occupation numbers ni| = A^, ni| = 0, 
not — 0; and tiqi = 0. Using Eq. (||) we see that the cor- 
responding fermion state has fermion occupaton numbers 
n-mi = 1 for TO = 1, • • • , A^^ and no| = , n„i| = 0, i.e. 
it is the fully polarized state with a single hole in the 
m — orbital. This state is the unique state in the one 
hole Hilbert speace with the maximum possible total bo- 
son angular momentum {M = A^) and the maximum 5*2 
(= A/2) The set of boson states with angular momentum 
M — N—SM and Sz = N/2~6Sz have boson occupation 
numbers satisfying 

noi + Tiot — SM 

noi+'nii = SS^. (4) 

For fixed 6M and SSz, the state may be specified by 
uqi, which can assume values from to the mimimum 
of 6M and SSz', the number of states is g{SSz,SA4) = 
1 + mi[5Sz, 6M]. We now deduce the total spin quantum 
numbers of the quasihole states from this expression. 

Since the Hamiltonian is spin-rotationally invariant, S'^ 
and Sz must be good quantum numbers, so that all eigen- 
states occur in spin-multiplets with degeneracy 25" -I- 1 for 
total spin S. Furthermore, assuming that edge effects are 
irrelevant, the Hamiltonian is invariant under simultane- 
ous translation of all coordinates. It follows that SM is 
also a good quantum number, and that each member of 
a spin-multiplet has associated with it a large orbital de- 
generacy which scales with the system size. Orbitally de- 
generate states can be generated from a seed state with 
minimum SM by repeated application of the operator 
which lowers the center-of-mass angular momentum [ p^ . 
For example, when Sz — N/2 this procedure generates 
holes which occur in successively larger single-particle an- 
gular momentum states. 

The Sz = N/2 — 1 manifold has one state with SM = 
and two states for each SM > 1. One state at each 
SM is the Sz = N/2 - 1 member of the S = N/2 spin- 
multiplet with the same SM. It follows that there is 
one S = N/2 — 1 spin-multiplet with an orbital degener- 
acy which is reduced by one compared to the S = A/2 
states. Continuing in the same way we may conclude 
that there is a single spin-multiplet with S = N/2 — K 
for each nonnegative integer K with orbital degeneracy 
N^ — K. Thus the quantized Skyrmion states occur in 
degenerate manifolds labeled by an integer K and with 
dimension (A'^^ — K){N — 2K -I- 1). The spin degeneracy 
(A — 2K + 1) is lifted by the Zeeman coupling and is 
the quantum counterpart of the arbitrary global orien- 
tation of a classical Skyrmion. The orbital degeneracy 
(A"^ — A') is the quantum counterpart of the arbitrary lo- 
cation of the center of a classical Skyrmion and is lifted 
by a disorder potential. The density of Skyrmion states 
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in the presence of Zeeman couphng and weak disorder is 
indicated schematically in Fig. ^ 

Repeated application of the total spin-lowering oper- 
ator and the center-of-mass lowering operator allows all 
states in the S = N/2 — K manifold to be generated 
from the seed state which has Sz ~ S and SM — K. 
To determine the many-body wavefunction of this spin- 
ier Skyrmion state we consider the effect of a total spin 
raising operator, 5'+, on the K + 1 boson states which 
occur at 5M = 5Sz = K: 

5*+ |noi = k, nof ~ K — k.nn = K ~ k) (5) 
1 

= E ^mT^™il"oi = k,noT = K - fc,nix ^ K ~ k) 

m=0 

« ^J{K - k)N\nai = k,no^ = K - fc,7iix = K - k-1). 

Here bl^ and &,„ are boson creation and annihilation op- 
erators. In the last form of Eq. (||) we have included 
only the m = 1 term which dominates for finite K and 
— > oo because ^7ii| « ^/N ^1. In this limit the state 
with k = K is annihilated by 5*+ and is therefore the seed 
state of the S = N/2 — K Skyrmion multiplet. In first 
quantization the (unnormalized) boson wavefunction of 
this state is 

E [ n U),exp(-|z,f/4)] 

il,---,iK ji£{iK} 

[ n |T);^;cxp(-|z,|V4)] (6) 

13{ik} 

In these wavefunctions, the sums are over distinct particle 
indices, up-spin quasiparticles occupy states with angular 
momentum m — 1 whereas down-spin particles occupy 
TO = states; this correlation between spin and angu- 
lar momentum states is an essential aspect of Skyrmion 
states [|,§. 

We now establish some relationships between the prop- 
erties of these wavefunctions and previously known re- 
sults. As explained above the K = Q Skyrmion wave- 
function is identical to a Hartree-Fock quasihole. On the 
other hand, the large K limit of the these wavefunctions 
may be related to classical field theory Skyrmions. Con- 
sider the fermion wavefunction 

oo 

|vI'(A)) = E A^l*^) 

i<j I 

exp(-|zHV4). (7) 

For large A the sum in Eq. (|^) will be dominated by terms 
in a relatively narrow range of K values |^(A)) 
is precisely the single-Slater-determinant proposed by 
Moon et al. M as a microscopic trial wavefunction for 



the Skyrmion on the grounds that, for large A, it gives a 
spin-texture in precise agreement with that of the classi- 
cal Skyrmion Q of size A. 

The more general Hartree-Fock (HF) single-Slater- 
determinant wavefunctions of Ref. |^ may be written in 
the form 

The HF approximation is equivalent to minimizing the 
energy of this wavefunction subject to a normalization 
constraint. An obvious deficiency of the Hartree-Fock ap- 
proximation is its failure to refiect known symmetries of 
the Hamiltonian. In particular, the Hartree-Fock wave- 
function is not an eigenstate of Sz or M . As discussed 
by Nayak and Wilczek this failure is readily reme- 
died by projecting the Hartree-Fock state onto a state of 
definite Sz (and therefore definite M). For the hard-core 
model, the HF equations may be solved analytically, 
with the result \um? = 1 - \vm? = >? /[X"^ + 2{m + 1)], 
where A is a free parameter. It is easily verified that 
the corresponding wavefunction is precisely |5'(A)). In 
the case of the hard-core model the projection of the 
Hartree-Fock wavefunction onto a state of definite Sz 
yields the exact Skyrmion wavefunctions, suggesting that 
this seemingly ad hoc procedure might be generically ac- 
curate. 

We remark that all the results we have obtained here 
to a hole in the ferromagnetic ground state apply equally 
well to the case of a particle added to the ferromagnetic 
ground state because of an exact particle-hole symme- 
try |]l8| which holds around v = \. For the case of the 
hard-core model the energy of the particle states is 4Vb, 
independent of K . Our analysis is also readily general- 
ized for fractional filling factors v = 1/to, although in 
that case the quasiparticle states cannot be generated 
by particle-hole transformation. For general interactions 
the Skyrmion energy is dependent on K] the minimum 
energy state may occur &t K — where Hartree-Fock 
theory is valid, aX K ^ oo where the classical field the- 
ory becomes valid, or at an intermediate value of K 
where explicit expressions for the energy are not avail- 
able. The many-particle wavefunctions which we have 
derived for the Skyrmion states are exact only for the 
hard-core model. The accuracy of these wavefunctions 
for general interactions, which has been confirmed by ex- 
act diagonaliztion jl^, rests on the existence of a gap 
for charged excitations and they are ,in this sense, analo- 
gous to Laughlin's trial wavefunctions for incompressible 
states at fractional filling factors. This work was sup- 
ported in part by NATO Collaborative Research Grant 
No. 930684, by the National Science Foundation under 
grants DMR-9416906, and DMR-9503814 and by CICyT 
of Spain under contract MAT 94-16906. HAF acknowl- 
edges the support of the A. P. Sloan Foundation and the 
Research Corporation. Helpful conversations with Steve 
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FIG. 1. Schematic single Skyrmion density of states. For 
quasihole states the orbital degeneracy for each S = N/2 — K, 
K = 0, 1,2, and Sz = —S,---,S is lifted by disorder 
producing a finite width band of states. The spin-multiplet 
structure persists in the presence of disorder. The energetic 
offset of bands with the same K and different Sz is due to 
Zeeman coupling. For the situation illustrated, the Zeeman 
spin-splitting energy is comparable to the disorder produced 
band width. The dependence of the Sz ~ S energy on K de- 
pends on the interaction Hamiltonian and the strength of the 
Zeeman coupling; the situation illustrated where the lowest 
energy state occurs at Sz = S = A'^/2 — 3, is typical. The 
Sz ~ S = N/2 — 5 band has been removed from this illustra- 
tion for clarity. 
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